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We present a theoretical investigation of the electronic structure of rutile (metallic) and Mi and 
M2 monoclinic (insulating) phases of VO2 employing a fully self-consistent combination of den¬ 
sity functional theory and embedded dynamical mean field theory calculations. We describe the 
electronic structure of the metallic and both insulating phases of VO2, and propose a distinct mech¬ 
anism for the gap opening. We show that Mott physics plays an essential role in all phases of VO2: 
undimerized vanadium atoms undergo classical Mott transition through local moment formation (in 
the M2 phase), while strong superexchange within V-dimers adds significant dynamic intersite corre¬ 
lations, which remove the singularity of self-energy for dimerized V-atoms. The resulting transition 
from rutile to dimerized Mi phase is adiabatically connected to Peierls-like transition, but is better 
characterized as the Mott transition in the presence of strong intersite exchange. As a consequence 
of Mott physics, the gap in the dimerized Mi phase is temperature dependent. The sole increase of 
electronic temperature collapses the gap, reminiscent of recent experiments. 

PACS numbers: 71.27.-|-a 71.30.-|-h 71.15.-m 74.20.Pq 


VO 2 is of great practical importance as its MIT is very 
close to room temperature, and its ultrafast switching 
between metallic and insulating phases can be used for 
device applications [l|. VO 2 is also of fundamental in¬ 
terest, due to the complex interplay between electronic 
correlations and structural distortions, and as a result 
the physical mechanism responsible for the gap forma¬ 
tion has remained under debate. 

Structurally, VO 2 undergoes a transition from a high 
temperature rutile phase 0 (R) to low temperature mon¬ 
oclinic phases. In the latter, under ambient pressure, the 
vanadium atoms dimerize and tilt with reject to the ru¬ 
tile c axis, giving rise to the Mi phase [3[. Due to the 
dimerization it is expected that the electronic states as¬ 
sociated with the overlapping d-orbitals along the rutile 
c axis (oig states) split into bonding-antibonding sub¬ 
bands. In addition, the antiferroelectric displacement 
leads to an upshift of the subbands associated with the 
remaining t 2 g states (ep. As a result, a gap would 
appear between the aig and subbands, suggesting 
thus that VO 2 undergoes a Peierls-type transition [^. 
However, density functional theory (DFT) calculations 
showed that these structural distortions alone cannot give 
rise to a gap in the Mi phase [1, @ . In addition, a dis¬ 
tinct monoclinic phase (M 2 ) appears when the system is 
under uniaxial stress or doped with Cr^+, Al^+, Fe^’*', or 
Ga^"*" [3,Q. In this phase half of the vanadium atoms pair 
along the rutile c axis, while the other half experience a 
zigzag-like distortion along the same axis. In particular, 
the existence of localized electrons in the zigzag chains 
also suggests that electronic correlations indeed play a 
role in the gap formation of VO 2 . 

The cluster-dynamical mean field theory (DMFT) has 
successfully described several aspects of the physics of 
VO 2 , but some discrepancies between different imple¬ 


mentations remain. Biermann et al. [9| showed that 
the Ml phase is insulating, with the gap of 0.6 eV in 
agreement with experimental findings [10l | . Furthermore, 
the authors found that electronic intersite correlations, 
within the vanadium dimers, renormalize down the aig 
bonding-antibonding splitting in comparison with DFT 
calculations. In their proposed mechanism the Mi phase 
can be viewed as a renormalized Peierls insulator. In con¬ 
trast, by means of ab initio linear scaling DFT-|-cluster- 
DMFT calculations, Weber et al. ll| observed that the 
gap formation of the Mi phase is mainly due to an 
orbital-selective Mott instability concerning the aig elec¬ 
tronic states in our notation. Looking at the occupancy 
of the 3d shell, they obtained at around two electrons 
per vanadium, resulting in four electrons per vanadium 
dimer, suggesting that the Mi phase is not a renormal¬ 
ized Peierls insulator Q. 

Important fundamental challenges remain: i) the close 
proximity of the Mi, M 2 , and R phase in the phase dia¬ 


gram ll2| calls for an unified framework in which these 


phases can be simultaneously described; ii) another fun¬ 
damental question is whether the MIT transitions will 
take place if we fix the structure and change only the 
temperature. 

It has not been possible to properly address these is¬ 
sues within an approach based on an effective model 
from which the oxygen degrees of freedom are elimi¬ 
nated, and which contains only V-3d electrons (Hubbard 
model). Indeed, earlier DMFT works on the Hubbard 
model [1, 13, 141 showed that for a given set of Hubbard 
U and J parameters, the gap in the Mi insulating phase 
is too robust, persisting to very high temperatures [^, 
while the mass renormalization in the R phase is too 
small compared to experiment [illil. Both of these 
effects result from placing this material too far from the 
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Mott transition boundary. Consequently, the M 2 phase 
was not described before with this approach, as half of 
V-atoms would not undergo Mott-Hubbard transition. 

In this letter we solved these problems by describ¬ 
ing VO 2 with modern all electron embedded DMFT ap¬ 
proach, where in addition to correlated vanadium atoms 
the itinerant states of oxygen are included in the Dyson 
self-consistent equation 15, 0- We present a compre¬ 
hensive picture and describe all the phases of VO 2 with 
the same values of the {U,J) parameters, including the 
M 2 phase which has not been considered previously in 
DMFT treatments. Mott physics is central for the proper 
description of all the phases even though the Mott insta¬ 
bility is arrested in the Mi phase. 

Our theory leads to a different physical picture for 
the gap opening in monoclinic phases of VO 2 , and most 
importantly the possibility of the collapse of the Mi 
insulating state by temperature. Our results indicate 
the presence of significant intersite correlations within 
the vanadium dimers, which in turn lower the aig sub¬ 
band in relation to the e^. In particular, we notice that 
nonlocal dynamic correlations enhance the aig bonding¬ 
antibonding splittiM, in contrast to what was reported 
by Biermann et al. [9|. The electrons in all phases of VO 2 
are in the near vicinity of the Mott transition, but the 
true pole in the self-energy, signalizing the local moment 
formation, occurs only in the paramagnetic M 2 phase on 
undimerized V-atoms and its aig orbital. In the Mi phase 
and in the antiferro-ordered M 2 phase, the singularity of 
the self-energy is arrested as the ordered states are adi- 
abatically connected to Peierls and Slater insulators, re¬ 
spectively. The adiabatic connection between the weakly 
and the corresponding strongly interacting states, makes 
the Mott mechanism hard to distinguish from alternative 
scenarios, nevertheless, collapsing a large insulating gap 
with raising electronic temperature is possible only in 
the Mott state in the presence of strong superexchange, 
as was found for example in the cluster-DMFT study of 
the 2D Hubbard model [l7|. Hence, according to our 
results the Mi phase is best characterized as the Mott 
phase in the presence of strong intersite superexchange 
within the V-dimers, while the undimerized V-atoms in 
the M 2 phase undergo canonical Mott transition asso¬ 
ciated with local moment formation in the presence of 
weak superexchange. 

Rutile and Mi phase. - We first address the R phase at 
temperature of 390 K. In Fig.jTfa) we show the calculated 
total, t 2 g and projected density of states (DOS). 

The R phase is metallic with strongly renormalized 
Og orbitals compared to corresponding DFT results (not 
shown). A weak lower Hubbard band (LHB) is observed 
at around -1.09 eV, and a broad upper Hubbard band 
(UHB) around 2.54 eV. The former is in good a gree ment 
with experimental photoemission measurements [l0| , and 
previous theoretical findings Q. The LHB and UHB are 
mainly of aig and character, respectively. 



FIG. 1 . (Color online) DFT-|-DMFT-based total (black 
dashed line) and projected DOS of (a) R and (b) Mi phase of 
VO2. The projections to aig, 6 ^( 1 ), eg( 2 ), and Cg states are 
shown in blue, red, green, and brown lines, respectively. For 
the Ml phase, the solid (dashed) bine line corresponds to the 
projection on the bonding (antibonding) aig molecular state. 


We next present the electronic structure of the Mi 
phase. Within the V-dimer, treated as cluster in DMFT, 
it is useful to adopt the symmetric and antisymmetric 
combination of orbitals. The associated self-energies are 
denoted as the bonding Tjb,a and antibonding Eat,a com¬ 
ponents, where a = {aig, 6^(1), 6^(2)} [l^. In Fig. [T{b) 
the total, t 2 g, and projected DOS are shown. At first, 
we notice the opening of a gap of 0.55 eV, between the 
aig and subband, which is in good agreement with ex¬ 
perimental and previous theoretical findings [§-ll|. As 
expected, upon dimerization and antiferroelectric distor¬ 
tion, the aig subband splits in bonding (solid blue) and 
antibonding (dashed blue) states, while the subband 
upshifts in comparison with the R phase. The aig bond¬ 
ing orbital has a coherent peak at around -0.30 eV, while 
the antibonding orbital has two incoherent peaks remi¬ 
niscent of LHB and UHB located at -1.5 eV and 2.58 eV, 
respectively. The coherent peak and the satellite below 
it were also observed in previous theoretical [i, [III and 
experimental 0 works. However, in our case, the UHB 
does not represent a weak coherent peak associated with 
the antibonding state, as in Ref. citebiermannPRL. The 
bonding-antibonding splitting (relative to the DFT and 
DFT -I- single-site DMFT calculations (not shown)) in¬ 
creases upon inclusion of intersite dynamic correlations 
within the dimer, in contrast to Ref. l9|. 


We also computed the optical conductivity of R, Mi, 
and M 2 phases of VO 2 , which are presented in Fig.[2{a). 
We calculate the plasma frequency of the R phase, de¬ 
fined as ujp = 2/(7reo) a(oj)doj, with A = 2eV and 
A = 0.1 eV, the latter denoted by At 390 K we ob¬ 
tained for the ratio u!p{LDA)/uip{DMFT) = 2.24, and 
ojp{DMFT) = 3.28 eV, in good agreement with Ref. [ll| . 
We notice a mild temperature dependence of this ratio, 
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FIG. 2. (Color online) (a) Real part of the optical conduc¬ 
tivity of R, Ml, and M2 phases of VO2. (b) Real part of 
intersite self-energies, on imaginary frequency axis, of aig-aig 
(black) and eg(l)-eg(l) (red) states of Mi phase, (c) Imag¬ 
inary part, on real frequency axis, of bonding (dashed lines) 
and antibonding (solid lines) self-energies associated with the 
aig (orange) and eg(l) (cyan) dimer electronic states. 


decreasing with decreasing temperature. This signals an 
increase of kinetic energy due to gain of coherence, char¬ 
acteristic of Mott-Hubbard system. On the other hand, 
the ratio oj*^{LDA)/oj*'^{DMFT) is much larger: at 
100 K it is 4.71, and increases as the Drude peak broadens 
(at 390 K it is 6.12). Finally, the DMFT mass enhance¬ 
ments, extracted from the slope of the self-energy, are 
m*Iruband = 3.5, 3, 2.5 for oig, e^(l), and e^(2), respec¬ 
tively. In the conventional Hubbard model picture, the 
mass enhancement {m*/rriband) and the ratios of plasma 
frequencies {ujp{LDA)/u!p{DMFT)) should be very simi¬ 
lar, and both are used to measure the correlation strength 
in the literature M- In VO 2 these measures of corre¬ 
lations differ because the plasma frequency ojp includes 
substantial contribution from the interband transitions, 
which strongly mix with less correlated oxygen and Cg 
states, resulting in apparently weaker correlations than 
given by m^ 2 g/rriband, a physics beyond Hubbard model. 
The low energy ratio of has negligible interband con¬ 
tribution, and at zero temperature it should be equal to 
mt 2 g/rriband, but at finite temperature it gets large due 
to broadening and the loss of strength of the Drude peak. 

Effects of nonlocal dynamic correlations. - Having de¬ 
scribed the electronic structure of both R and Mi phases 
of VO 2 , we next investigate how nonlocal dynamic cor¬ 
relations contribute to the opening of a gap between oig 
and Cg subbands in the Mi phase. Fig. [^Dc) shows the 
bonding and antibonding components of the imaginary 
part of the self-energy associated with the aig and eg(l) 
dimer electronic states. 

We notice that once the Peierls instability occurs in 
our calculation, the Mott instability is arrested, hence 


there is no pole in the imaginary part of the self-energy 
associated with the aig or Cg states. This excludes an 
orbital-selective Mott-Hubbard mechanism as the source 
of the gap of the Mi phase, as proposed in Ref. [ll|. Our 
results bear strong resemblance with the Mott transition 
of the Hubbard model in cluster-DMFT [^, where the 
local singlet state of the cluster dominates the low energy 
properties of the model. The energy gain to form the 
strong bonding state on the cluster is here not just due 
to increased hopping between the two V-atoms, but it is 
primarily due to the gain of the exchange energy, which 
is stronger than kinetic energy, as the latter is strongly 
reduced due to proximity to the Mott transition. The 
precise calculation of this exchange energy is beyond this 
work since it requires the computation of the momentum 
dependent spin-susceptibility which in turn depends on 
the two particle vertex function. 

Next we define the local and intersite self-energies, 
which are related to the self-energies in the dimer basis as 
j^iocai ^ + and = i(Ef,-E,b) 0. 

In Fig. HKb) we also show the real part of the inter¬ 
site components of self-energies associated with the same 
electronic states as those in the inset. We observe that 
the component of eg(l) is negligible, but, notably, one 
can see that depends strongly on the fre¬ 

quency in the low-energy part. This indicates the pres¬ 
ence of strong intersite electronic correlations within the 
vanadium dimers, which in turn lower the aig bond¬ 
ing state. As a result, the bonding-antibonding split¬ 
ting increases and a gap between the aig and Cg ap¬ 
pears (for more details about the bonding-antibonding 
splitting see the supplemental material). Thus, one can 
consider that nonlocal correlations give rise to an effec¬ 
tive oig — oig frequency dependent hopping ta^g-aig + 
ReT,’ff^_g^^^{ioj), which properly takes into account the 
aig bonding-antibonding splitting. A similar observation 
was previously proposed for the low-temperature phase 
of Ti 203 [2l|, but required a strong intersite Coulomb 
interaction for opening the gap. 

Metallization due to hot carriers.- More recently, nu¬ 
merous experimental studies have reported the existence 
of monoclinic-like metallic phases of VO 2 [ii-il. This 
transition was also induced by the application of fem¬ 
tosecond laser pulses on VO 2 films. As pointed out 
by Wegkamp et al. [ 2 ^, the photoexcitation gives rise 
to hot carriers, which have an associated temperature 
much higher than the lattice temperature. Motivated by 
this fact, we performed calculations considering a much 
higher temperature, i.e. T = 900 K, for electrons in the 
Ml phase. 

In Fig. [SKa) and (b) we show the calculated spectral 
function and the projected DOS of Mi at 332 and 900 K, 
respectively. From the spectral function at 900 K, one 
can see the closing of the gap. In particular, we ob¬ 
serve that aig and Cg (1) subbands shift towards the Fermi 
level, the latter shifting more than the aig subband. To 
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FIG. 3. (Color online) Spectral function and projected DOS 
of Ml phase at 332 (a) and 900 K (b). In (c) the spectral 
function of M 2 phase at 332 K is shown. 


understand why the gap vanishes at 900 K, we exam¬ 
ine the temperature dependence of the self-energies. In 
Fig.SKa) we present the real part of intersite for 

these two temperatures. In addition, in the inset we show 
the antibonding component of Re^ab,aig and the bond¬ 
ing component of i?eE{, on the real frequency axis. 
Here, we observe that in the low-energy 

part, and ReTiab.aig are strongly suppressed with increas¬ 
ing temperature. Therefore at 900 K the renormalization 
of the aig subband decreases significantly in comparison 
to that at 332 K. Thus, this subband is shifted towards 
the Fermi level. In relation to the 6^(1) subband, we note 
an enhancement of the ReT,b,e^(i) from -0.46 eV to 0.77 
eV, leading to an downshift of this subband. 

M 2 phase.- Our cluster-DMFT simulation of the M 2 
phase, using the same values of U and J employed for 
the other phases, are presented in Fig. [3Kc). The size of 
the insulating gap is 0.58 eV. In addition we show in 
Fig.[2Ka) the calculated optical conductivity of M 2 phase. 
Interestingly, we observe that both monoclinic phases are 
very similar, providing further support for the Mott point 
of view. 

For this simulation we allowed anti-ferromagnetic or¬ 
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FIG. 4. (Color online) (a) (ita) at 332 K 

(dashed lines) and 900 K (solid lines) of Mi phase. Inset: 
Re'Eab,aig{bj) (blue) and Re'£ii,e^ (red) at 332 (dashed 
lines) and 900 K (solid lines), (b) ReT,aig{u}) of undimer¬ 
ized V-atoms in the antiferromagnetic phase of M 2 phase, (c) 
of bonding and antibonding a\g self-energies of V- 
dimers (orange lines) and a\g states of undimerized V-atoms 
(black lines). 


dering of undimerized vanadium atoms, but we find a 
similar gap also in paramagnetic simulation of this phase. 
In both simulations, we find that the aig orbital on the 
dimerized V-atoms is renormalized due to intersite corre¬ 
lations, in similar way as in the Mi phase. However, the 
splitting mechanism of the aig orbital on undimerized V- 
atoms depends on the magnetic ordering. In particular, 
in the anti-ferromagnetic phase the aig orbital is renor¬ 
malized by the spin-dependent real part of the local aig 
self-energy, which depends strongly on frequency in the 
low-energy part as shown in Fig. HKb). In contrast, in 
the paramagnetic phase the aig self-energy has a pole in 
the imaginary part as can be seen in Fig. [UJc). This in¬ 
dicates that in the paramagnetic phase, the aig orbital 
of the undimerized V-atoms splits due to the canonical 
Mott instability, proving that the M 2 phase should be 
characterized as a Mott insulator. In addition, this Mott 
instability also leads to an orbital polarization in favor of 
the aig orbital of undimerized V atoms by about 0.34/V 
atom with respect to the same orbital in R phase. It is 
noteworthy that the Cg orbital occupation decreases al¬ 
most by the same amount (-0.28/V atom). Finally, we 
mention that this orbital polarization leads to the stabi¬ 
lization of the aig electronic states and therefore coop¬ 
erates for the stabilization of the insulating state in the 
M 2 phase. 

Conclusions.- Our work highlights the importance of 
describing simultaneously all phases of a strongly cor- 
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related material. The importance of Mott physics in all 
the phases of VO 2 was stressed early on in the pioneering 
work of Pouget and Rice 27|. This physics now emerges 
from a quantitative first principles method, and its impli¬ 
cations for many physical quantities has been elucidated. 
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METHOD 


We performed all-electron charge self-consistent 
DFT-fDMFT calculations, with implementation de¬ 
scribed in Ref. [^. This implementation avoids down¬ 
folding, and rather extremises DMFT functional in real 
space. This implementation has been tested in nu¬ 
merous materials, such as transition metal oxides 
iron-based superconductors [3-0, actinides [lol - fl^ , lan¬ 
thanides [H-lil, and transition metal compounds [ibl . 


17l| . Here it is applied, for the first time, to a system 


that requires cluster-DMFT treatment. 


In our calculations we included dynamic correlations 
within single-site and also within two-site cluster on 
dimerized V-atoms (in monoclinic phases). In partic¬ 
ular, the single-site DMFT is very accurate for non- 
dimerized vanadium atoms (see Refs. Bui) , while 
the dimerized sites require cluster-DMFT treatment. 
Within our calculations, we obtain the electronic charge 
density, chemical potential, impurity levels, self-energy 
(S), and Green’s function self-consistently. We men¬ 
tion that the full charge self-consistency considerably im¬ 
proves on one-shot calculation, as has been for exam¬ 
ple shown in the case of metal-insulator transition in 
V 2 O 3 [ 2 ^. The all-electron DFT calculations were per¬ 
formed within Perdew-Burke-Ernzerhof generalized gra¬ 
dient approximation (PBE-GGA) [^, as implemented 
in Wien2K code [ 2 ^. The quantum impurity problem 
was solved by means of Gontinuous time quantum Monte 
Garlo (CTQMG) calculations [ 2 ^, at fixed distinct tem¬ 
peratures for R and monoclinic phases. To compute the 
self-energy on the real frequency axis we performed an¬ 
alytic continuation using the maximum entropy method. 
Throughout our calculations we employed the Coulomb 
interaction U = 6.0 eV and Hund’s coupling J = 1.0 eV. 
We have also used the local axis system proposed by Ey- 
ert 1 ^ . wherein the d,,; 2 _y 2 orbitals mediate the u-type 
d — d overlap along the rutile c axis, while the remain¬ 
ing t 2 g orbitals, i. e. dxz and dyz, mediate the 7r-type 
d — d overlap. In our notation, the former is denoted as 
aig and the latter as the doubly degenerated states. 
In particular, in our investigation {6^(1), 6^(2)} denotes 
the {dxz, dyz} and {dyz, dxz} orbitals in R and Mi phase, 
respectively. In our calculations only the t 2 g states were 
considered in the impurity model, while the remaining 
Bg states were treated within DFT (GGA). Finally, to 


investigate the R and Mi phase we used the experimen¬ 
tal structural parameters from Ref. [i^ and Ref. [23 |. 
respectively. 


DYSON EQUATION OF THE SOLID 

The correlations are quite local in real space, but they 
are not very local in a Wannier representation, which 
eliminated oxygen degrees of freedom and builds Hub¬ 
bard model, keeping vanadium d-bands only. The DMFT 
(c-DMFT) approximation, which treats exactly the dy¬ 
namic correlations on a give atom (or cluster), is not 
a very good approximation in the latter case, because 
Wannier functions are quite delocalized and extending on 
oxygen sites. On the other hand, DMFT is an excellent 
approximation in real space, when dynamic correlations 
are expanded in terms of quasi-localized atomic orbitals, 
i.e., 

E,^(r,r')= X (r|d>';)($(;|E|$C;,)(<I>(;,|r') (1) 

mm' 

where |d)((j) form a complete basis set around vana¬ 
dium atoms centered at yt, and m, m' enumerate atomic 
degrees of freedom. Our quasi-atomic orbitals van¬ 
ish outside the radius Rmt ~ The single¬ 

site DMFT approximation truncates the self-energy to 
(d>((,|S|$C;,) = (d>((,|S|<I>^,), and cluster-DMFT 

keeps intersite terms ($fl, |S|$(1,) within the given clus¬ 
ter (cellular-DMFT (2^). As the functions |d)((j) are very 
localized, the negligence of off-site terms on un-dimerized 
vanadiums is an excellent approximation. 

When dynamic correlations are expanded in the form 
of Eq. [U one can not build a low energy (Hubbard-like) 
model containing only a small number of bands, because 
all bands hybridize with atomic degrees of freedom. In¬ 
deed, the Dyson equation 

Giu:{r,r') = {{iuj -f /r -f - I4-s(r))(5(r - r') 

-S,4r,r'))“' (2) 

couples all bands with the atomic-degrees of freedom, 
and therefore the Dyson equation has to be solved either 
in real space, or, equivalently in complete basis, such as 
Kohn-Sham basis. This approach is also called embed¬ 
ded DMET, as the self-energy is embedded into the large 
Hilbert space, rather than the more common truncation 
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of the degrees of freedom of the solid to a Hubbard model 
through downfolding. 

The truncation to local Feynman diagrams (which is 
the DMFT approximation) is less severe approximation 
in such formulation of the theory than in the alternative 
Hubbard model formulation. Another advantage of such 
formulation of DFT+DMFT is that the approximation is 
conserving and stationary [ 2 ^ . as it is derivable from a 
Luttinger-Ward (LW) functional. Namely, LW functional 

r[G] = Trlog(-G) - Tr((Go-i - G-^)G) + 

Eh+Xc[p] + ^[Glocal] — EH+Xc[Plocal\ (3) 


is stationary at the DFT+DMFT solution, i.e., 
(5r[G]/(5G = 0 at G = Gdmft- 
To prove that, we take the derivative of the above func¬ 
tional, and obtain 


_ 1^-1 _ (^-1 
SG~ ° 


Vks 


(5$[G, 


local\ 


SG 


(4) 


where we denoted ^[Giocai] = ^[Giocai] - En+xclpiocai]- 
This derivative vanishes if the self-energy in the Dyson 
equation (Eq. [5|) satisfies 


^ _ SGiocal 5^[Giocal] - . 

5G ~ 6G SGiocai ^ ’ 

In our embedded formulation of DMFT, Giocai is a pro¬ 
jection of the Green’s function to the set of quasi-atomic 
orbitals 

GioU^.v') = ($(;,|G|$C;,) i^CW). (6) 

mm' 

/i, G cluster 


Using this form of the projector, we can compute 
5Giocai/SG, and hence the form of the self-energy at the 
stationary DFT+DMFT solution: 

E(r,r') = ^(r|ci>(^) (<,|r')(7) 

OUlocal 

mm 

fi, fi' G cluster 

This is exactly the form postulated in Eq. [U hence we 
proved that DET+DMET set of equations are derivable 
from the F functional, and are stationary. Note that 
within DMFT approximation, is computed by 

solving an appropriate quantum impurity model. Note 
also that in deriving Eq. [7] we assumed that the basis 
functions do not depend on the Green’s function 

(or the charge density), and hence /5G = 0. In 

the alternative downfolding implementations of charge 
self-consistent DET+DMET, the functions are re¬ 
placed by the Wannier functions, which essentially de¬ 
pend on the charge density of the solid, and there¬ 
fore on the Green’s function. In such Wannier imple¬ 
mentations, the derivative of the projector (originating 
in nonzero 5 I'hm) /5G) has always been neglected, and 
therefore charge-self consistent DET+DMET method is 
non-stationary in Wannier implementations. 


LOCAL AND INTERSITE SELF ENERGIES 

Within our approach the cluster self-energy is diago¬ 
nal in the momentum basis, with quantities Ek=o and 
E,k=tv The cluster self-energy in the site representation, 
on the other hand, contains the local and intersite 
self-energies in the following form 


Se = 




Using the cluster symmetry one can notice that En = 
S 22 and Ei 2 = E 21 , where the former is the local self¬ 
energy and the later the intersite self-energy. By means 
of a Fourier transform one can obtain the self-energy in 
the momentum basis from the local and intersite self¬ 
energies [ 2 ^ . 









( 8 ) 


where Nc is the number of sites within the cluster, K is 
the cluster momentum, and Rij denote the sites within 
the cluster. In our cluster we have R — 0,1 and K = 0, tt. 
Thus, one obtain that 


E_fs:=o = Ell + Ei2, (9) 

and 

U_fs:=7r = Ell — U 12 . (10) 


Based on the molecular scenario with bonding and an¬ 
tibonding orbitals, we change the names of and 

E,k=tv to bonding (Ef,) and antibonding (E^f,) self¬ 
energies, respectively. Therefore, the two previous equa¬ 
tions lead to the following linear relations: 

Ell = i(E6 + E,6), (11) 


Ei 2 = i(E6-E,6), (12) 

which are used in the main part of the text. 

As we showed in the main text, the imaginary part 
of the self-energies associated with the bonding and an¬ 
tibonding states in the Mi phase vanishes at the Fermi 
level. Hence, we can discuss the bonding-antibonding 
splitting within an effective band structure picture, in 
which the effects of electronic correlations are taken into 
account expanding the self-energies to first order. In this 
case the bonding and antibonding bands are renormal¬ 
ized by the quantities 





gjREfc 

duj 




-1 


(13) 
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and 


Zab 



duj 


uj—O/ 


-1 


(14) 


The value of these quantities are well defined since there 
is no pole at zero frequency in the Matsubara self-energy. 
As a result, we can obtain the energies of renormalized 
bands 


Eb = Zb{eb + ^Eb{uj = 0)), (15) 


Eab — Zabi^ab + 3?Sab(w — 0)), (16) 

where et and Cab are the energies of noninteracting bond¬ 
ing and antibonding bands, respectively. 

Therefore, the bonding-antibonding splitting energy, 
A = Eab — Eb, can be written as 


A =Zab{eab + 3fiSafa(0)) - Zbiet + 3fiEfa(0)) 

= '^- - -^(Cab + 3iSab(0) + eb + 5RE;,(0)) 

+ {Zab + Zb) (^^^ + ^^^^(0) _ _ 3^^^(0))^ (^7^ 


Taking the bonding and antibonding self-energies in 
terms of Sn and E 12 we obtain 


Zb 


_ QjRSii I x , / ag?Si2 
1-*- dui 1 duj 

(-1 _ dWSii I \2 _ r dWSi2 
dui 1 dui 



(18) 


Zab = 


A da; (. doj 

('I _ QjRSii I '12 _ /' 0j)?Si2 
1-*- doj l‘^=0l I duj 



(19) 


Notice that since Zab, Zb, and (1 — U=o) are pos¬ 
itive, it is the sign of U=o that determines which 

orbital is more correlated. 

Finally, on the basis of equations (fT51) and (ITOl) the 
bonding-antibonding splitting energy is then given by 


A =- 


( 1 - 


doj 


j=o) 


( 1 - 


OSRSii 


.= 0 )^ - ( 


asEi 


duj V 0 ,^ 

X (eab-eb-2SRSi2(0)) 

r ausi2 1 '1 

V Oa, k=0j 


.= 0 )' 


raSREi 


Cl dMl .11 I „'|2 _ Ci_ 

1-*- doj l‘^=0l V dui 

X {eab + Cb + 25RI]ii(0)). 


0=0^ 


( 20 ) 


This formula has a simple physical interpretation. The 
first term is just the “local” Ziocai multiplied by the bare 
splitting, enhanced by the intersite self-energy, while the 
second term arises from the asymmetry of the local crys¬ 
tal field levels corrected by local self-energy, and would 
vanish in the particle-hole symmetric limit. 
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